Encoded Recoupling and Decoupling: An Alternative to Quantum Error Correcting 
Codes, Applied to Trapped Ion Quantum Computation 
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A recently developed theory for eliminating decoherence and design constraints in quantum com- 
puters, "encoded recoupling and decoupling", is shown to be fully compatible with a promising 
proposal for an architecture enabling scalable ion-trap quantum computation [D. Kielpinski et al., 
Nature 417, 709 (2002)]. Logical qubits are encoded into pairs of ions. Logic gates are implemented 
using the S0rensen-M0lmer (SM) scheme applied to pairs of ions at a time. The encoding offers 
continuous protection against collective dephasing. Decoupling pulses, that are also implemented 
using the SM scheme directly to the encoded qubits, are capable of further reducing various other 
sources of qubit decoherence, such as due to differential dephasing and due to decohered vibrational 
modes. The feasibility of using the relatively slow SM pulses in a decoupling scheme quenching the 
latter source of decoherence follows from the observed 1/f spectrum of the vibrational bath. 

PACS numbers: 03.67.Lx,89.70.+c,42.50.Lc 



INTRODUCTION 



In the quest to construct a scalable quantum computer, 
a recent proposal Q advocating an array-based approach 
to quantum computing (QC) with trapped ions appears 
particularly promising. Ions are stored for later process- 
ing in a memory region, then transported to an inter- 
action region, where pairs are coupled in order to enact 
quantum lògic gates. This proposal overcomes some of 
the design constraints that plagued the original Cirac- 
Zoller (CZ) ion-trap-QC proposal 0], which prevented 
the latter from becoming fully scalable. In the new pro- 
posal an encoding of a single logical qubit into the states 
{|0l) = I-IT), = ITJ-)} 01 t wo trapped- ion (physical) 
qubits is used. Quantum lògic gates are implemented us- 
ing the S0rensen-M0lmer (SM) scheme p, |J (see also 
related schemes by Milburn et al. || |^), which has 
the advantage of reduced sensitivity to motional state 
heating compared to the CZ proposal. The encoding 
into {| J,|), | tl)} is useful because these states form a 
decoherence-free subspace (DFS) [[?J [|, ^, [ï(| with rc- 
spect to collective dephasing, a process whereby the envi- 
ronment introduces identical random phase modulations 
on groups of physical qubits |Tïj . In the context of the 
"quantum charge coupled device" (QCCD) proposed in 
, such a process is one of two dominant sources of de- 
coherence. The DFS encoding reduces the collective de- 
phasing problem by several orders of magnitude ]Ï2| . A 
method to perform universal QC using the SM scheme on 
these DFS qubits was proposed in jïj , and independently 
in ||. 

The DFS encoding {| U)} is well known |, 

|Ï4| , ft5[ pi , and its utility against collective dephasing 
was dcmonstrated experimentally using photons |l7| and 
trapped ions [ jl8| . The notion of universal QC using a 
DFS has been explored in general in {[(], [Ï9], M 
for Hamiltonians that always preserve the DFS, in p2, 
B3L p4| using a combination of DFS and active quan- 



tum error correction methods |25|, ||, |7|, gg |(| |ï), 
and in |3^] , using an approach wherein transitions out of 
the code space are suppressed by continuous observation. 
Still more generally, the notion of universal QC while 
overcoming decoherence as well as design constraints has 
been explored by us and coworkers at a theoretical level in 



a series of recent papers 
This theory uses a combination of qubit encoding into 
a DFS with selective recoupling |34| and dynamical de- 




coupling 42, M 



|5C| ] operations, 
so we refer to it as "encoded recoupling and decoupling" 
(ERD). The utility of ERD as a general method for quan- 
tum simulation, universal QC, and decoherence suppres- 
sion has also been stressed and explored by Viola [|TJ . 

The main focus of ERD so far was on sòlid state 
ü M |35|, |36|, |37l |38l |§ g ^ and NMR §| QC 



proposals. Here we present a unified treatment of the 
ERD ideas, and show that they apply also in an atòmic 
physics setting, namely the QCCD ion-trap proposal [jïj. 
Specifically, we show here how to perform universal QC 
on DFS qubits comprised of pairs of trapped ions, by 
using the SM scheme for quantum lògic, in a manner 
that involves manipulating only pairs of ions at any given 
time, while always perfectly preserving the DFS encoding 
(see Brown et al. J53[ for an interesting alternative set of 
ion pairs-only lògic gates, which, however, does not pre- 
serve the DFS at all times). By applying strong and fast 
dynamical decoupling ("bang-bang" |Ï2|, Q) SM pulses 
we show how to further drastically reduce sources of de- 
coherence beyond collective dephasing. While a qubit 
is being stored an example of such a source of deco- 
herence is deviations from the collective dephasing ap- 
proximation. While a qubit is being manipulated for the 
purposes of QC, coupling to vibrational modes is nec- 
essary [||, ||, ||, ||, |6|, and decoherence of these vibra- 
tional modes due to patch-potential noise is the second 
dominant source of qubit decoherence |LE 
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A method to suppress vibrational mode decoherence (as 
well as heating, which is not as serious a problem when 
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the SM scheme is used), employing a version of the dy- 
namical decoupling method known as "parity kicks" , was 
proposed and discussed in detail by Vitali & Tombesi 
(VT) [[Ï6| This method uses a fast+strong modula- 
tion of the trapping potential. We present here an alter- 
native decoupling method for suppressing decoherence of 
ion-trap qubits due to their coupling to decohered vibra- 
tional modes, that operates directly on the qubit (spin-) 
states. The feasibility of this scheme, in spite of the rela- 
tive slowness of the SM pulses, follows from the observed 
1// spectrum of the vibrational bath |5^, p7] . The con- 
centration of most of the bath spectral density in the 
vicinity of the low, rather than the high-frequcncy cut- 
off, implies much relaxed constraints on the decoupling 
pulses compared to those usually assumed JÏ2], ^(j, [I?) . 
A full analysis of this result will be presented in a forth- 
coming publication f5g| . 



More generally, we show here how all sources of deco- 
herence beyond collective dephasing can in principle be 
suppressed using sufficiently strong+fast SM pulses. This 
includes so-called bath-induced "leakage errors" , wherein 
the system-bath coupling induces transitions into or out 
of the qubit subspace |37], [59|, We provide feasi- 

bility estimates for the decoupling pulses and find that 
they are within current experimental reach. The overall 
picture emerging from this work is that ERD provides 
a means for a robust, decoherence-resistant implementa- 
tion of universal QC with trapped ions. Experimental 
implcmentation of the ERD method should be possible 
with current ion-trap technology and we suggest a few 
experiments. 

The structure of the paper is as follows. In section |ü] 
we review the DFS encoding into two spins and the asso- 
ciated lògic gates. We show how our previous formulation 
thereof can be reinterpreted in the context of acting on 
pairs of trapped ions within the SM scheme. Wc also 
present a method for coupling pairs of encoded qubits 
using pulses that involve controlling only pairs of ions at 
a time, while always preserving the DFS encoding. In 
the subsequent sections we discuss how to reduce deco- 
herence. In section III we review the decoupling method, 
emphasizing its application to trapped ion arrays. We 
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then proceed to apply the ERD method: in section 
we show how to eliminate the residual differential dephas- 
ing contribution to decoherence using SM pulses; and, in 
section [v] we discuss how to reduce all further sources of 
decoherence, including the component that arises due to 
coupling to decohered motional states. Then, in section 
VI we show how to fully implement ERD, i.e., we show 
how to combine universal QC via recoupling over DFS- 
encoded qubits with decoherence suppression via encoded 
decoupling. To makc ERD fully effective for trapped 
ions we suggest to combine it with the VT potential- 
modulation method. Concluding remarks are presented 



II. ENCODED UNIVERSAL LOGIC GATES IN 
ION TRAPS 



To fix tcrminology we first connect the methods devel- 
oped in p3l |34|, [55| to the gates proposed for trapped ions 
in [Q . Let Xi, Yi , denote the Standard Pauli matrices 
af,af,a^, acting on the ith physical qubit (we will use 
both notations interchangeably) . In 3^] it was shown 
that for the code {\0 L ) = ||T), |U) = fi)} the encoded 
logical operations (involving the first two physical qubits) 
are 



~{X 1 X 2 +Y l Y 2 ) 



X12 

Y\2 = — \YxX2 — X{Yí 
Z12 = 77(^1 — Z 2 ). 



(D 



These operations form an su(2) àlgebra (i.e., we think 
of them as Hamiltonians rather than unitary operators). 
We use a bar to denote logical operations on the en- 
coded qubits. In |33|, |3^] these logical operations 
were denoted by T" 2l a £ {2, y, z}, and a detailed analy- 
sis was given on how to use typical solid-state Hamilto- 
nians (Heisenberg, XXZ, and XY models) to implement 
quantum lògic operations using this DFS encoding. E.g., 
the term X\X 2 + Y{Y 2 is the spin-spin interaction in the 
XY model, and Z\ — Z 2 represents a Zeeman splitting. 
A static Zeeman splitting and a controllable XY interac- 
tion can be used to generate a universal set of lògic gates, 
a result that has very recently been applied in the con- 
text of spin-based QC using semiconductor quantum dots 
and cavity QED |}9j . Similar conclusions hold when the 
XY interaction is replaced by a Heisenberg |54| |6C], [nj or 
XXZ interaction |35| , or even by a Heisenberg interaction 
that includes an anistropic spin-orbit term |Ï8). We re- 
mark that, as first shown in ]R| [HI , the various types of 
exchange interactions can be made universal also without 
any single-qubit terms (such as a Zeeman splitting), by 
encoding into three or more qubits |H| |6^, |6^, |64|, p 3 k a 
result that has been termed "encoded universality" Kjq|. 



A. Logic gates on two ions encoding a single logical 

qubit 

S0rensen and M0lmer proposed a quantum lògic gate 
that couples two ions via a two photon process that vir- 
tually populates the excited motional states of the ions 
H] . The SM scheme works well even for ions in thermal 
motion, while the CZ scheme requires cooling the ions to 
their motional ground statc. The SM scheme involves ap- 
plying two làsers of opposite detuning 5 to the two ions. 
Ideally the Lamb-Dicke limit should be satisfied: 



in section VII. 



(n+ í)i] 2 < 1, 



(2) 
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where rj is the Lamb-Dicke parameter and n is the mean 
vibrational number. Deviations from the Lamb-Dicke 
limit lead to fidelity reduction that is proportional to rf 
|^| . The time required to prepare a maximally entangled 
statc using the SM scheme is 



? 7 



K 



(3) 



where Í7 is the Rabi frequency and K is an integer [Q. 
For realistic parameters, in the strong field limit (K = 1 
in Eq. (12) of |Q), tsm can be made as short as 1/xsec. 

In JjJ it was shown that the SM two-ion gate can be 
expressed as follows. The unitary gate 1/2(0, <f>%, </> 2 ) was 
introduced, which we here rename Uij(9, <pi, <j> ): 



Uij(9,4>i,4>j) = expfyOXfaXfy) 



cos Oljj + i sin BXfaXç 



(4) 



where 



Xé = X cos < 



Y sin < 



The phase <f>i is the phase of the driving làser at the zth 
ion, while 9 oc Í2 and can be set over a wide range of 
vàlues (J, ^7). Introducing the operators 

X l3 = h&iXi - YYj), Y 3 = ±(Y t Xj + XiY-) (5) 



(denoted Rf 3 , R?- respectively in [|33j, |34|, |35|) we can 
express 



Uij(0,<f>i,(j)j) = cos 61 + i sin (cos A<f>ijXij 

+ sin A(j)ij Yij + cos í>y Xij + sin Y i3 ) , 

(6) 



where $ 



i+<pj ■ It is simple to check that X i3 and Y i3 
annihilate the code subspace {|0l) = |4T)>|1l) = líi)} 
and have non-trivial action (as encoded X and Y) on 
the orthogonal subspace {|||), |TT)}- Therefore, as also 
observed in |Q and Q , upon restriction to the DFS we 
can write: 

f/y (Mi, Üij(9, A^ij) 

= ex.p(i9X A4>tj ) = cos 61 + i sin 9X A 4 >lj ■ 

(7) 

The fact that Uij depends only on the relative phase A(j>ij 
is crucial: this quantity can be controlled by adjusting 
the angle between the driving làser and the interatomic 
axis, as well as by small adjustments of the trap voltages 
(which, in turn, control the trap oscillation frequency, 
and hence the ion spacing), whereas it is much hardcr 
to control the absolute phase & Jï|, |6^, |6^], and hence 
also This is why the code subspace {| J.T), | |J.)} 

enjoys a preferred status over the subspace {|I4), |TT)}- 
A thorough theoretical analysis of the approximations 



leading to the gate Uij(0,<fii,(f)j) is given in [|| (see also 
|Ï3f for an abbreviated exposition that emphasizes the 
connection to computation in a DFS). 

Let us establish the connection between the seemingly 
distinct sets of lògic operations in Eqs. (|ï|),(|^). To do 
so, we only need to use Eqs. @,(@) while ignoring the 
component that annihilates the DFS (X, Y). Then: 

exp(i0X 12 ) = U 12 (6,4>,<t>) = Üi2 (6,0) 
exp(i0F 12 ) = [/ 12 (0,^0+|) = Í7 12 (0,|) 

exp(i6>Zi 2 ) = exp(i^y"i 2 ) exp(zSXi 2 ) exp(-ï^Fi2) 
= C/ 12 ( J , tt/2) C7 12 0) É7 12 ( - J , tt/2) . (8) 

The third linc follows from the elementary operator iden- 
tity 

X = X cos 4> + Y sin (j> = e -^z/2 Xel· <pz/2 ^ 

which holds for any su(2) angular momentum set 
{X,Y,Z}, i.e., operators that satisfy the commutation 
relation [X, Y] = 2iZ (and cyclic permutations thereof), 
in particular also the encoded operators {X, Y, Z}. 

This proves the equivalence of the two sets of operators. 
Using these results and Eq. (^), a more direct connection 
can be written in terms of the Hamiltonians: 



X12 
Y 12 



X^X,), 

X<j,X ( j )+7r /2i 



(10) 

(11) 



where the equivalence is meant in terms of a projection 
of the RHS Hamiltonians to the DFS. In the context of 
ion-trap QC the lògic gate Ü(9, A0) can be performed 
directly, so it may be more intuitively useful than the 
{X, Y, Z} set. Eq. (|J) shows that by properly choosing 
9 and A^íj all single DFS-encoded qubit gates can be 
performed. 



B. Entangling gate between pairs of encoded 
qubits involving four ions 

In Q the following unitary gate was introduced 



Ui = exp^i-X^X^X^X^) 

= (h h hh-iX^ X^ X 4>3 X <t>i ] 
^ -/| (^2^34 - iX A4>i2 X A<j>34 ) 
= exp(-ZjX A(/ , i2 XA0 34 ). 



(12) 



Th is g ate, also considered in slightly less general form 
in Q, can be used to entangle two DFS-qubits. It in- 
volves simultaneous control over two phase differences 
A<f> 12 , A034, and thus control over the motion of two pairs 
of ions. The case A(j) 12 = A0 3 4 = was used in |6é| 
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to demonstrate entanglement of four trapped-ion qubits, 
but this choice is not unique. 

We now come to an important point that was not ad- 
dressed in Q: in order for the DFS encoding to offer 
protection against collective dephasing during the exec- 
tion of the entangling gate, collective dephasing condi- 
tions must prevail over all four ions. To see this, notc 
that a differential dephasing term such as [Z\ — Z3) <g) B 
(where B is a bath operator) does not commute with t/4, 
so that if such a term exists during gate execution then 
the DFS will not be preserved, according to a theorem 
in pjj . On the other hand, collective dephasing over 
all four ions, expressed by a system-bath coupling of the 
form Q2í—i Zi) ® B, does commute with U4, so that in 
this case the DFS is preserved |Ï9[] . While deviations 
from collective dephasing over pairs of ions have been 
shown experimentally to be small ||Ï8| , this may not be 
the case over the length seales involving four ions [ jOj . 
We discuss in section [V how to create such extended 



collective dephasing conditions. 

Taken together, the results in this section show how 
universal QC can be implemented using trapped ions by 
applying the SM scheme to pairs of ions at a time, each 
encoding a DFS qubit. The DFS encoding takes care 
of protecting the encoded information against collective 
dephasing. We now move on to a discussion of how to 
reduce additional source of decohcrcncc. 



III. D YN AMIC AL DECOUPLING PULSES AND 
THEIR APPLICATION TO TRAPPED IONS 

Let us briefly review the decoupling technique, as it 
pertains to our problem (for an overview see, e.g ., [ pï l ). 
Decoupling, as proposed by Viola and Lloyd |f|[ [43|, rc- 
lies on the ability to apply strong and fast pulses, in a 
manner which effectively averages the system-bath inter- 
action Hamiltonian Hsb to zero. A quantitative anal- 
ysis was first performed in for pure dephas- 
ing in the linear spin-boson model (which corresponds 
to the ohmic case of the Caldeira-Leggett model [|69[): 
Hsb — ® B, where B is a Hermitian boson oper- 
ator. The analysis was recently extended to the non- 
linear spin-boson model, with similar conclusions about 
performan ce |49|| . Imperfections in the pulses were con- 
sidered in ]4q], and it was shown that an optimal value 
for the pulse period can be found. Since the decou- 
pling pulses are strong one ignores the evolution under 
Hsb while the pulses are on, and since the pulses are 
fast one ignores the evolution of the bath under its free 
Hamiltonian Hb during the pulse cycle. The simplest 
example of eliminating an undesired unitary evolution 
U = expf— ü(Hsb + Hb)], is the "parity-kick sequence" 
p2| , [IM p6[ . Suppose we have at our disposal a fully 
controllable interaction generating a gate R such that 
lí R conjugates U": RÏUR = W . Then the sequence 
UR^UR = I serves to eliminate U . A simple example of a 
parity kick sequence is the following. Assume we can turn 



on the single-qubit Hamiltonian ÇÏXj for a time 7r/2í7. 
This generates the single-qubit gate Xj — i exp(— i^Xj). 

Suppose that H SB = EÜi T, a e{x, y ,z} 7?°"? ® B ï ■ Each 
term in Hsb either commutes or anti-commutes with Xj . 
If a term A in Hsb anticommutes with Xj then the 
evolution under it will be conjugated by the gate Xj-. 
Xjey^{-iAM)X,j = cxp(-iXjAXjAt) = exp(iAAt). 
This allows for selectively removing this term using the 
parity-kick cycle, which we write as: [Aí, Xj, At, Xj]. 
Reading from right to left, this notation means: apply 
Xj pulse, free evolution for time Aí, repeat. Suppose 
that we can also apply the single qubit gate Yj. Then, 
since every system factor in the above Hsb contains a 
single-qubit operator, it follows that we can selectively 
keep or remove each term in Hsb by using the parity- 
kick cycle. Note, however, that without additional sym- 
metry assumptions, this procedure, if used to eliminate 
all errors, requires a number of pulses that is exponential 
in the number of qubits N The reason is that 

without symmetry assumptions we will need at least two 
non-commuting single-qubit operators per qubit (e.g., X , 
Y) , and we will need to concatenate decoupling pulse se- 
quences. Below we show how to dynamically generate 
such symmetries, in a way that avoids this exponential 
sealing (for a discussion of this point see the Conclusions 
section). 

Note that in the analysis of the parity kick cycle we 
ignored Hsb and Hb while R was operating; this is jus- 
tified by the strength assumption. The bath Hamiltonian 
Hb commutes with the applied pulses, but its effect is 
very important since [Hb,Hsb] 7^ in general. Thcrc- 
fore if the bath has spectral components at freqüències 
higher than the inverse of the interval between decou- 
pling pulses, then the bath density matrix will pick up 
phases that are essentially random, and this effect will 
show up as decoherence (for a quantitative analysis see 
j42], £Ï7|, |8[ H| . Hence it is commonly assumed that 
the pulse interval, Aí, should be small compared to the 
inverse of the high-frequeney cutoff lü c of the bath spec- 
tral density I(ui) [ïj|, which also sets the scale of the 
bath-induced noise correlation time í c (the speed assump- 
tion). It can be shown that the overall system-bath cou- 
pling strength "fsB is then renormalized by a factor Atu) c 
after a cycle of decoupling pulses ]Í4| , or that the bath- 
induced error rate is reduced by a factor proportional to 
(Aí/í c ) 2 Q. Using a Magnus expansion ]70fl , it can be 
shown that there is a hierarchy of decoupling schemes, 
whereby jsb is renormalized by a factor (Aíw c ) fe , where 
fe > 1 is the order of the decoupling scheme [H . The im- 
plication for single-qubit dephasing, Hsb = ^1sbZ ® B 
(B is a dimensionless bath operator), is that the dephas- 
ing time T2 is increased by a factor 1 / ( Aíwc) 2 * [fn] . Thus 
it seems crucial to be able to apply pulses at intervals 
Aí <C 1/uj c . However, as shown first by Viola & Lloyd 
p2| , and then by VT in their quantitative analysis of 
a vibrational mode linearly coupled to a boson bath, a 
finite-temperature bath sets another, thermal timescale 
that must be beat in order for the decoupling method to 
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work 



Let the system-bath coupling be 



H 



SB 



(13) 



where a (b) is an annihilation operator for the system 
(fcth bath) vibrational mode, and 7 is the (for simplicity 
uniform) energy damping rate. In the context of trapped 
ions the bath is provided by fluctuating patch-potentials 
(due, e.g., to randomly oriented domains at the surface 
of the elèctrodes or adsorbed materials on the elèctrodes) 
p6[ . Then VT showed that the decoupling pulse interval 
(in fact, the entire cycle time) must be shorter also than 
the thermal decoherence time 

ídcc(T) = { 7 [l + 2n(T)]}-\ 

where n(T) = [ e huJ o/k B T _ -q-i j g mean vibrational 
number of the system oscillator at thermal equilibrium 
with temperature T, and lüq is the frequency of the os- 
cillator, i.e., the system is described by the harmònic os- 
cillator Hamiltonian Hs = hw^a^a. Thus the timescale 
condition for successful decoupling is 

Aí<min{l/w c ,í dcc (T)}. 

As shown in the VT analysis, the timescale ídec(T) is es- 
pecially relevant for vibrational mode decoherence in ion 
traps, which as already mentioned above, is responsible 
for qubit decoherence during quantum lògic gate opera- 
tions. 

However, for trapped ions experimental evidence so far 
points to a 1//" spectrum for the vibrational bath over 
a range 1 — 100MHz J56|, p.5], implying that there is no 
clear high-frequency cutoff lo c . Encouragingly, in a re- 
cent experiment involving a charge qubit in a small su- 
perconducting elèctrode (Cooper-pair box), a version of 
parity-kick decoupling was successfully used to suppress 
low-frequency energy-level fluctuations (causing dephas- 
ing) due to í/f charge noise This suggests that 

decoupling can help even in the absence of a clear cutoff 
frequency. Recent theoretical results support this ob- 
servation |{3é|: for í/f noise most of the bath spectral 
density is conccntratcd in the low, rather than the 
high-end of the frequency range. 

In spite of the apparent 1/ f a spectrum in trapped ions, 
VT used a cutoff estimate of lü c < 100 MHz @, and 
showed that suppression of vibrational decoherence can 
be accomplished by pulsing the oscillation frequency lüq 
of the ion chain (i.e., by pulsing the trapping potential), 
provided Aí < í/lu c ~ lnsec, a nd T < lOmK. 

Given the estimate in Section II A of t$m ^ for the 
SM gate, it is clear that we cannot hope to satisfy the 
strict Aí < lnsec timescale requirement which would be 
needed in order to use decoupling directly on the qubit, 
rather than the vibrational modes, assuming the VT es- 
timate of lü c . However, the theoretical analysis |Q and 
the success of parity-kick decoupling in the presence of 
1/f noise in the charge qubit case [Í72| suggests that it 



may well be worthwhile to apply decoupling pulses on 
the qubit in addition to, or perhaps instead of, the VT 
trapping-potential-modulation scheme. 

Now let us comment on the strength assumption. Here 
we must make sure that the amplitude of the decoupling 
pulses, íl, can be made much stronger than the system- 
bath interaction 7 in Eq. (|T^), i.e., the heating rate from 
the vibrational ground state of the ion chain. Experimen- 
tal measurements of 7 are very sensitive to trap geomctry, 
secular frequency, and si ze J pq| , and range from a few Hz 
to a few tens of KHz [[Ï6| |73||~On the other hand, one can 
have fi as high as 1MHz Wjï, so the strength assumption 
can be comfortably satisfied. This does come at a price, 
however, since in the strong field limit the SM gate is per- 
turbed by a term that yields direct, off-resonant coupling 
of the qubit | f ) and | J.) states without changes in the 
vibrational motion j|. This is a unitary gate error that 
decreases the gate fidelity by (N / 2) (fi / S) 2 , where N is 
the number of ions participating in the gate Table II] . 
This forces us to be in a parameter regime where 57 <C S. 
In principle it is possible to exactly cancel this effect if 
the duration of the làser pulses is chosen so that both 
Eq. (H) and the condition tsm = K'n/S are satisfied, 
where K' is an integer and S is the detuning. However, 
in the context of decoupling we will also need to satisfy 
conditions such as íirgM = 7r/m where m is an integer. 
Putting these conditions together yields 



K'ir 



n/m 



S = mK'tt 



Çl — \pK — Tr/m =>• n = msÍK 
77ÍÍ 

While there is no problem with the first of these, the 
second condition implies that we cannot be in the Lamb- 
Dicke limit, Eq. (||). Therefore exact cancellation is not 
possible in our case, and we must resort to Q <C S in 
order to keep the fidelity reduction to a minimum. On 
the other hand, the kind of unitary error that is caused by 
off-resonant coupling can be corrected by optimal control 
pulse shaping methods ffq |. 

Finally, we note that fluctuations of various experi- 
mental parameters, such as intensity and phase fluctu- 
ations of the exciting làsers, can cause pure dephasing 
of the vibrational modes, in addition to the dissipative 
coupling described above J7(| . Clearly, the success of de- 
coupling strategies hinges on strong suppression of such 
fluctuations, as in the threshold theorem of fault tolerant 
quantum error correction p5| , p6| , |30| , pTf . 

To conclude, the discussion in this section indicates 
that the experimental viability of decoupling schemes in 
ion traps is rather promising, although it is hard to es- 
timate their success at this point. In the following sec- 
tions the analysis will be carried out at a more abstract 
level, emphasizing the algebraic conditions for a success- 
ful implementation of ERD. In the end it will be up to 
an experiment to decide the usefulness of the proposed 
schemes. 
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IV. CREATING COLLECTIVE DEPHASING 
CONDITIONS USING DECOUPLING PULSES: 
REDUCING DECOHERENCE DURING 
STORAGE 

One of the important advantages of the DFS encoding 
{Hi), |T!)} is that it is immune to collective dephasing. 
However, other sources of decoherence inevitably remain. 
In this and the following section, we algebraically classify 
all additional decoherence effects and show how they can 
be eliminated, in particular in the context of trapped 



A. Creating collective dephasing on a pair of ions 

First, let us analyze the effect of breaking the collec- 
tive dephasing symmetry, by considering a system-bath 
interaction of the form 



or, in ion-trap terms: 



fJ -dcph(2) v 



-iHsi 



r £/i2(-^,0)e-^ r [/ 12 Q,0) 



-í(Z 1 + Z 2 )®B* o1 2t 



Bf 



Bo 



2' ' L ^2' 

(14) 

where [Zy (9 , A0y) was defined in Eq. (R), and we used 
the identification found in Eq. (||). This equation means 
that the system-bath coupling effectively looks like col- 
lective dephasing at the end of the pulse sequence. Thus, 
the system is periodically (every 2r) projected into the 
DFS. 

In order for the the procedure described in Eq. ( fÏ4| ) 
to work, the SM gate Í7i2(±|,0) must be executed at 
a timescale faster than the cutoff frequeney associated 
with the fluctuating magnètic fields causing the differen- 
tial dephasing term in iï s . c ^ ph ' 2 ' ) . This cutoff has not yet 
been characterized experimentally, but the decay rate of 
the DFS-encoded state of two ions has been measured to 
be 2.2KHz [Q. Using this as a rough estimate for the 
cutoff frequeney, we see that the procedure of Eq. ( |Ï4| ) is 
likely to be attainable with fast (tsm ~ 1m s ) SM pulses. 



where Bf,B 2 are arbitrary bath operators. This de- 
scribes a general dephasing interaction on two qubits, 
and we can expect this to be the case during storage of 
trapped ion qubits in the QCCD proposal. The source of 
such dephasing during storage is long wavelength, ran- 
domly fluctuating ambient magnètic fields p8|, that ran- 
domly shift the relative phase between the qubit | |) and 
| |) states through the Zeeman effect. The interaction 
can be rewritten as a sum over a collective dephasing 
term Z\ + Z 2 and another, differential dephasing term 
Z\ — Z 2 , that is responsible for errors on the DFS: 



H -dcph(2) _ 
n SB 



= (Z x + Z 2 ) ® B z col + (Z 1 - Z 2 ) ® B z àii . 



Here B* ol = (Bf + Bf) /2 and Bfa = (Bf -Bf)/2. If 
Sfjj were zero then there would only be collective dephas- 
ing and the DFS encoding would offer perfect protection 
j77[ . However, in general B^ if ^ 0, and the DFS encoding 
will not suffice to offer complete protection. 

The crucial observation is that, since Z± — Z 2 oc Z\ 2 
[recali Eq. ( pi)], the offending term causes logical errors 
on the DFS ]37j]. Then the problem of B* i{ ^ can be 

solved using a series of pulses that symmetrize ' 
such that only the collective term remains, as shown in 
p(i|, p>H (7§j. To do so note that since the offending term 
cx Z12, it anticommutes with X X2 = 3(^1^2 + ^1^2) ■ 
At the same time X12 commutes with Z\ + Z 2 - This al- 
lows us to flip the sign of the offending term by using a 
pair of ±7r/2 pulses in X12, while leaving only the col- 
lective term. Evolution with the flipped sign followed by 
unaltered evolution leads to cancellation of the offending 
term. Specifically |36|: 



e- lHsBT e~ 



íHsbt 



-i{Z 1 +Z 2 )®Bl al 2r 



B. Creating collective dephasing on a block of four 



So far we have discussed creation of collective dephas- 
ing conditions on a sin glc DFS qubit. However, as men- 
it is essential for the reliablc exe- 



tioned in Section II B 



cution of an entangling lògic gate to have collective de- 
phasing over all four ions participating in the gate, even 
if only two are coupled at a time. A procedure for creat- 
ing collective decoherence conditions over blocks of 3, 4, 6 
and 8 qubits was given in Here we show how to do 
the same for a block of 4 qubits with collective dephasing. 

Let us start with a general dephasing Hamiltonian on 
iV ions, and rewrite it in terms of nearest-neighbor sums 
and differences: 



H 



deph 
SB 



N 

^Z t ®B t 

i=l 
N/2 

( z 2j + Z 2j -i) ® B+ + (Z 2j - Z 2j -i) ® B-j, 



where Bfj 



(B 2 j ± _B 2 j-i)/2. As noted above, Z : 



oc Z 2 ,_i 



2j'-l,2j, 



2j - 

so that to eliminate all nearest- 



neighbor differences of the form (Z 2 j 
the collective decoupling pulse X nn ■ 



Z 2 j-i) we can use 



\x-> 



p ~iH SB r x iH SB r x -\ _ p -i2TZ , yí(Z 2j +Z 2j - 1 )®B+ 

or, in ion-trap terms: 
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N/2 



3 = 1 



Ü 2j - 1>2j {-~,Q) 



-íHsbt 



N/2 



3=1 



>U 2j - l<2j {~,0) 



The next step is to eliminate next-nearest neighbor differential terms. To this end let us rewrite the outcome of the 
X nn pulsc in terms of sums and differences over blocks of four ions: 

N/2 N/2 

( Z 2j + Z 2j -i) ® B+ = J2 [Z2 3 +2 + Z 2j+ i + Z 2j + Z 2j -i] <8 B+- + 

3=1 3=1 

N/2 

+ [(^3+2 - Z 2j ) + (Z 2j+1 - Z 2j _ x )\ ® B+f, 

3=1 

where SÍ' = (-B^ +2 ± j3^)/2. The term in the first line contains only the desired block-collective dephasing over 4 
ions. The term in the second line contains undesired differential dephasing terms that we wish to eliminate. But these 
terms once again have the appearance of encoded Z operators, between next-nearest neighbor ion pairs. Therefore 
we need to apply a second collective pulse X nnn = ^fj^ e l ^ X2j ~ 1 ' 2j+1 e l ^ X2: '· 2j+2 , that applies encoded X operators 
on these ion pairs. At this point we are left just with collective dephasing terms on blocks of 4 ions, as required: 



-i2T Y. "ií (Z2Í+Z23-! )®b + 



N/2 



3=1 



l u 2j ^ 2j+í (- - , 0) U 2j , 2j+2 (- 1 , 0) 



)®B. 



N/2 



3 = 1 



' U 2j - lt2j+1 (|, 0)U 2j}2j+2 {^, 0) 



(15) 



This pulse sequence is important to ensure that collective 
dephasing conditions will prevail during the execution of 
lògic gates between DFS qubits, as emphasized in Section 
IIB . 

To conclude, the procedures discussed in this section 
provide a means for engineering collective dephasing con- 
ditions in an ion trap experiment. We propose here to 
implemcnt these symmctrization schemes experimentally. 
Moreover, we propose to combine them with the lògic 
gates described in Sec. |j How to do this efficiently is 
discussed in Sec. VI below. 



V. REDUCTION OF ALL REMAINING 
DECOHERENCE ON A SINGLE DFS QUBIT 
DURING LOGIC GATE EXECUTION 

The reduction of differential dephasing errors, as in the 
previous subsection, is particularly relevant for storage 
errors. However, this is only the first step. Additional 
sources of decoherence may take place during storage, 
and in particular during the execution of lògic gates, the 
most dominant of which is qubit decoherence due to cou- 
pling to decohered vibrational modes, as discussed above. 
It is useful to provide a complete algebraic classification 
of the possible decoherence processes. This will allow us 



to see what can be done using SM-decoupling pulses. To 
this end let us now write the system-bath Hamiltonian 
on two physical qubits in the general form 



H 



SB - 

where 

-Hdfs = Span{ 



Hi 



Lcak 



h 



Logi 



H 



DFS 



HLeak 
HLogi 

Z 



ZI + IZ XY + YX XX -YY 
2 ' 2 ' 2 ' 
ZZ, II} 

Span{X7, IX, YI, IY, XZ, ZX, YZ, ZY} 

XX + YY YX-XY 
Span{X = ,Y = , 



ZI-IZ 



} 



(16) 



where I is the identity operator, XZ = X\Z 2 (etc), and 
where Span means a linear combination of these opera- 
tors tensored with bath operators. The 16 operators in 
Eq. ( |Ï6| ) form a complete basis for all 2-qubit operators. 
This classification, first introduced in |37| , has the fol- 
lowing significance. The operators in Hdfs either vanish 
on the DFS, or are proportional to identity on it. In ei- 
ther case their effect is to generate an overall phase on the 
DFS, so they can be safely ignored from now on. The op- 
erators in HLoak are the leakage errors: terms that cause 
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transitions between states inside and outside of the DFS. 
A universal and efficient decoupling method for elimi- 
nating such errors, for arbitrary numbers of (encoded) 
qubits was given in [^9| . Finally, the operators in PLogi 
have the form of lògic gates on the DFS. However, these 
are undesired lògic operations, since they are coupled to 
the bath, and thus cause decoherence. 

It is worthwhile to already emphasize that the oper- 
ator YI + IY G iÏLcak is of particular importance: As 
shown in Eq.43], this is the operator that describes 
qubit decoherence due to motional decoherence during 
application of the SM gate. 

In the previous subsection we showed how to eliminate 
the logical error Z, but we see now that this was only 
one error in a much larger set. To deal with the addi- 
tional errors it is useful at this point to introduce a more 
compact notation for the pulse sequences. We denote by 
[r] a period of evolution under the free Hamiltonian, i.e., 
U(t) = cxp(—iHsB T ) = [t], and further denote 

P = tfi 2 (-|,0) = expH|x 12 ). 

Thus Eq. ( |Ï4| ) can be written as: 

ex P H(5í + Bi)(Zi + Z 2 )t] = [r, P, r, pt]. 

Notice that this is an example of a "parity-kick" pulse 
sequence. 

As a first step in dealing with the additional errors, 
note that the symmetrization procedure [r, P, r, P^] can 
in fact achieve more than just the elimination of the dif- 
ferential dephasing Z\ — Z 2 terra. Since X\2 also anti- 
commutes with Yvx — 5(^1-^2 — -X1Y2) € PLogi, if such 
a term appears in the system-bath interaction it too will 
be eliminated using the same procedure. 

So far we have used a 5 X12 pulse. Interestingly, the 
Hamiltonian X±2 can also be used to eliminate all leak- 
age errors [|?]. To see this, note that J7i 2 (±7r, 0) = 
exp(±Í7rAi2) = Z\Z 2 - This operator anticommutes with 
all terms in PLeak- Hence it too can be used in a parity- 
kick pulse sequence, that will eliminate all the leakage 
errors. In particular, this pulse sequence will eliminate 
qubit decoherence due to motional decoherence, i.e., the 
error YI + IY e P Lc ak- 



At this point we are left with just a single error: 
2 ® B itself, in PLogi- Clearly, we cannot use a pulse 
generated by X12 to eliminate this error. Instead, to 
deal with this error we need to introduce one more pulse 
pair that anticommutes with A12, e.g., exp(±í-fYi2) = 

t/l2(±f,f). 

Let us now see how to combine all the decoherence 
elimination pulses into one efficient sequence. First we 
introduce the abbreviations 



II = Pi 2 (±7T,0) = eX P (±Í7TAl2 

Q = Ü 12 { 



IÚ = PP 



7T 7T 7T — 

_,-)=ex P H-Fi2) 



A = Pi 2 (±7t, -) = exp(±z^yi 2 ) =At =QQ (17) 
As argued above, the 7r pulse II eliminates PLeak^ 

exp[-z(P Logi + P DFS )2r] = [r, II, r, IT]. 

This may be sufficient in practice, since as argued above 
this pulse sequence eliminates the YI + IY term, and the 
DFS encoding takes care of collective dephasing. Thus 
we expect that using cycles of two pulses we can almost 
entirely eliminate the two most important sources of de- 
coherence. This expectation of course depends on the 
time scale requirement for decoupling being satisfied, as 
discussed in detail in Section III above. In practice it 



may well be advantageous to combine the DFS encoding 
and [t, II, r, II] pulse sequence with the VT method of 
pulsing the trapping potential JÏ6|, . 

Now let us discuss adding the extra pulses needed to 
achieve full decoherence elimination. The ir/2 pulse P 
eliminates Y and Z in PLogi- Combining this with the 
sequence for leakage elimination we have the sequence of 
4 pulses: 



-i(JWf*®B*)4r _ [u( t )TIU{t)IL]P^[U(t] 
= [r,n,r,P,r,n,r,pt], 



np(r)n]p 

(18) 



(where we have used IIP^ = P, IIP = PT). 

If we wish to entirely eliminate decoherence then we 
are left just with getting rid of the logical error due to 
X. To eliminate it we now combine with the y-direction, 
ir/2 pulse, Q: 



J 



[U(r)IlU(r)PU(r)UU(T)P 1 }Q' í [U(T)ILU(T)PU(T)IlU(r)P i: ]Q 
[t, n, r, P, r, II, r, P\Q\t, II, r, P, r, II, r, P\Q] 



(19) 



r 



which takes ten pulses. Unfortunately it is not possible 
to compress this further, since P^Q = (iX)(—ïY) = iZ 
and PtQT = — iZ, neither of which cannot be gen- 
erated directly (in one step) from the available gate 



Üij(8, A(f>ij) = cos 61 + i sin 6X a^- • Finally, note that 
in principle the last pulse sequence is applicable also to 
other QC proposals, such as NMR and quantum dots. 

One important caveat (mentioned in Section [Fn| above) 
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is that, because we need very strong and fast pulses, our 
gate operation may become imperfect. Specifically, off- 
resonant coupling and deviations from the Lamb-Dicke 
approximation may become important. The former in- 
troduces a term XI + IX into the Hamiltonian generat- 
ing the Uij(9, (f>i,(f>j) gate Sec. IIIA], which can cause 
unitary leakage errors from the DFS. These can in turn 
be reduced using the methods in 75, |7Í|. Whether the 
decoupling method we have proposed offers an improve- 
ment will have to be put to an experimental test. 



VI. COMBINING LOGIC GATES WITH 
DECOUPLING PULSES 

So far we have discussed computation using the en- 
coded recoupling method (Section , and encoded de- 
coupling (Sections |y|,|v|). We now put the two together 
in order to obtain the full ERD picture. At least two 
methods are available for combining quantum comput- 
ing operations with the sequences of decoupling pulses 
we have presented above. For a general analysis of this 
issue see 



A. Fast + Strong Gates Method 

The decoupling pulse sequences given in Sec. fv] "stro- 
boscopically" create collective dephasing conditions at 
the conclusion of each cycle. As noted above, this is 
equivalent to a periòdic projection into the DFS. This 
property allows for "stroboscopic" quantum computa- 
tion at the corresponding projection times ]2Ï| . Here 
the computation pulses need to be synchronized with 
the decoupling pulses, and inserted at the end of each 
cycle. The amount of time available for implcmcntation 
of a lògic gate is no more than the bath correlation time 
t c = 2it/lú c . Assuming the dominant decoherence contri- 



butions not accounted for by the DFS encoding to come 
from differential dephasing (sctting the r c time-scale) and 
1/ f noise, and that we already assumed that we can use 
pulses with interval Aí <C t c , it is consistent to assume 
that we can then also perform lògic gates on the same 
time scale. 



B. Fast + Weak Gates Method 

There may be an advantage to using fast but weak 
pulses for the lògic gates, while preserving the fast + 
strong property of the decoupling pulses. To sec how to 
combine lògic gates with decoupling in this case, let us 
denote by Hs = X^X^. the controllable system Hamil- 
tonian that generates the entangling gate Uij(8,(j)i,4>j) 
[recali Eq. (Q)]. Suppose first that we turn on this logic- 
gate generating Hamiltonian in a manner that is ncithcr 
very strong nor very fast, so that the system-bath interac- 
tion is not negligible while Hs is on (this obviously puts 
less severe demands on experimental implcmcntation). 
Then the corresponding unitary operator describing the 
dynamics of system plus bath is: 

Ü(t) = exp[-it(H 8 + H SB + H B )]. 

Now, if we choose Hs so that it commutes with the de- 
coupling pulses, then we can show that after decoupling 



Ü(t) ^ exp[-i2t(H s + H B )}, 



(20) 



provided t is sufficiently small. Tracing out the bath 
then leaves a purely unitary, decoherence-free evolution 
on the system. To prové this, assume we have cho- 
scn t' and the decoupling Hamiltonian H' s so that (i) 
exp(— it'H' s )HsB exp(it'H' s ) = —Hsb (the parity kick 
transformation) , and (ii) [H' s ,Hg] — 0. Then 



U{t)e- lt ' H 'sU{t)e 



= Ü(t)e 



-it[íís+e _ " H SHsBe u H S+H B ] 



e -it(H s +HsB+H B ) e -it(Hs-H S B + H B ) 
p -{2it{Hs + H B )+t 2 ([H S B,Hs] + [HsB,H B ])+0(t :s )} 



r 



where we have used the Baker-Campbell-Hausdorff for- 
mula, exp(aA) exp(aB) = exp{a(A + B) + ^-[A,B] + 
0(a 3 )}. 

Setting Hs = QS and Hsb = IsbS 1 (8) B we have 
the condition t -C in order to be able to 

neglect the 0(t 2 ) term [H S b,H s ]. Using ü - 1MHz, 
7ss ~ lOKHz as in Section [H, we fmd t <C 



10/isec. 

However, the more stringent constraint comes from the 
[Hsb, Hg] term, since Hb is not bounded for a harmònic 



oscillator. A more careful analysis then shows the famil- 
iar conclusion, that the bath should not be allowed to 
evolve for longer than its correlation time |^2|, [47| . 
Hence the actual requirement may still the far more strin- 
gent condition t <C í/ui c < lnsec for the decoupling pulse 
interval; see Sec. III. This cannot be satisfied with SM 
pulses, but in this case we can resort to the VT poten- 
tial modulation method. When we do this in conjunction 
with SM decoupling pulses we can be sure that Eq. (|20|) 
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is an excel·lent approximation. On the other hand, the 
rcquircments for al// bath spectral density are far less 
stringent and may be satisfied even with SM pulses alone 
|pH . Furthermore, for the rotation angle 9 = Çlt describ- 
ing the computation we have 9 yJVL/^sB < 10, which 
means that there is no restriction on applying large ro- 
tations. 

Let us now show how to efhciently combine lògic op- 
crations and decoupling pulses. For simplicity considcr 
only the case where we can neglect the X error, i.e., our 



decoupling sequence is the 4-pulse one given in Eq. (18). 
Suppose we wish to implement a logical X operation, i.e., 
exp(— i9Xi2). Recali [Eq. (|Ï0|)] that this involves turn- 

ing on the Hamiltonian Hg = ÇlxX^X^ I S S Q,xX\2 
between two physical qubits. Because the decoupling 
pulses P = exp(— i^Xi 2 ) and II = exp(±Í7rATi 2 ) are gen- 
erated in terms of the same Hamiltonian, they commute 
with Hg while eliminating Hsb (except for the terms 
in Hsb that have trivial action on the DFS). Thus the 
conditions under which Eq. ( p0| ) were shown to hold are 
satisfied. This allows us to insert the lògic gates into the 
four free evolution periods involved in the pulse sequence 
of Eq. Thus, the full pulse sequence that combines 

creation of collective dephasing conditions with execution 
of the lògic gate is: 



e -it(QxX l2 +H DFS ) _ ü(t/4)UÜ(t/4:)PÜ(t/4:)nÜ(t/i)P^ , 

(21) 

with Ü(t) = exp[-it(iT^ + H S b + H B )], and which, 
using the DFS encoding, is equivalent to the desired 
exp(— i9X\ 2 ). This involves 8 control pulses, 4 of which 
are of the fast+strong type (those involving P and II), 
and 4 of which must be fast, but need not be so strong 
that we can neglect Hsb- 

If we wish to implement logical Y operation, i.e., 
cxp(— i6Y 12), then we cannot now use P and II, since 
they anticommutc with Y yi and will eliminate it. In- 
stead we should use decoupling pulses generated in terms 
of Y12, which will also have the desired effect of elimi- 
nating iÏLcak, as well as X and Z logical errors, while 
commuting with the Y lògic operations (and for this rea- 
son can of course not eliminate Y errors). These are just 
the Q and A pulses defined in Eq. (|Ï7|). In ion trap terms 
this implies [recali Eq. (|ïl"|)] turning on the Hamiltonian 

Hg — ÇlyX^X^+K j2 1— > ÇlyY Yi between two physical 
qubits. Thus: 



e -ií(Ovri2+HDP S ) = Ü{t/A)KÜ{t/A)QÜ{t/A)MJ{t/A)Q\ 

(22) 

with U(t) = exp[—it(H^ + Hsb + Hb)\, and which, 
using the DFS encoding, is equivalent to the desired 
exp(-i0Fi 2 ). 

Finally, to generate single DFS-qubit rotations about 
an arbitrary axis we can combine Eqs. (|2l"|) , (p2|) ac- 
cording to the Euler angles construction. Given that 
Eqs. (j2Ï|),(p2|) each take 8 pulses, the Euler angle method 
will generate an arbitrary DFS-qubit rotation in at most 
24 pulses. 



Concerning gates that entangle two DFS-qubits, the 
situation is more involved, since now the next-nearest 
neighbor pulses in Eq. (^5|) , that create the collective de- 
phasing conditions on four ions, do not all commute with 
the Ui gate of Eq. (|Ï2]). Therefore here we must resort 
to the strong + fast method of the previous subsection, 
i.e., we need to synchronize the U4 pulses with the end 
of the decoupling pulse sequence. 

Taken together, the methods described in this section 
provide an explicit way to implement universal QC using 
trapped ions in a manner that offers protection against 
all sources of qubit decoherence, using a fast + strong 
(or fast + weak) version of the SM scheme, possibly in 
combination with the VT potential modulation method. 



VII. DISCUSSION AND CONCLUSIONS 

We have proposed a method of encoded recoupling and 
decoupling (ERD) for performing decoherence-protected 
quantum computation in ion traps. Our method com- 
bines the S0rensen-M0lmer (SM) scheme for quantum 
lògic gates with an encoding into ion-pair decoherence- 
free subspaces (each pair yielding one encoded qubit), 
and sequences of recoupling and decoupling pulses. The 
qubit encoding protects against collective dephasing pro- 
cesses, while the decoupling pulses symmetrize all other 
sources of decoherence into a collective dephasing inter- 
action. The recoupling pulses are used to implement en- 
coded quantum lògic gates, either during or in between 
the decoupling pulses. All pulses are generated directly 
using the SM scheme. We have provided numerical esti- 
mates of the feasibility of our scheme, which seem quite 
favorable. In order to achieve full protection against all 
decoherence it may be necessary to supplement ERD 
with the potential modulation method due to Vitali & 
Tombesi, in order to reduce vibrational mode decoher- 
ence. However, it may be worthwhile to test ERD with- 
out potential modulation first, as a significant reduction 
in decoherence can already be expected according to the 
results presented here. This is so because the vibrational 
bath has been found experimentally to have a l// a spec- 
tral density j5?j] , and there exists evidence that in such a 
case decoupling may be possible under moderate timing 
contraints 1531. 



As mentioned in Section [II, the dynamical decou- 
pling method requires an exponential number of pulses 
if the most general form of decoherence is to be sup- 
pressed, that can couple arbitrary numbers of qubits to 
the environment (total decoherence [^0[). This exponen- 
tial sealing is avoided here by focusing on decoherence 
climination inside blocks of finite size (e.g., at most four 
ions) where arbitrary decoherence is allowed. However, 
we have implicitly assumed that there are no decoher- 
ence processes coupling different blocks. This is a rea- 
sonable assumption for trapped ions, where the different 
blocks can be kept sufficiently far apart until they need to 
be brought together in order to execute inter-block lògic 
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gates. When this happens, ERD can still be efficiently 
applied on the temporarily larger block. 

It may be questioned whether there is any advantage 
in using ERD compared to methods of active quantum 
error correcting codes (QECC). Both ERD and QECC 
are capable of dealing with arbitrary decoherence pro- 
cesses, and are fully compatible with universal quantum 
computation. There are two main advantages to ERD: 
First, we need only two ions per qubit, compared to a 
redundancy of five ions per qubit to handle all singlc- 
qubit errors in QECC |27j]. So far experiments involving 
trapped ions have used up to four ions floq] , so that this 
encoding economy is a distinct advantage for near-term 
experiments. Second, our method is directly compatible 
with the SM scheme for lògic gates in ion traps. On the 
other hand it is not clear how to directly use SM gates 
for QECC. These are general features of ERD: economy 
of encoding redundancy and use of only the most easily 
controllablc interactions. On the other, the disadvantage 
of ERD compared to QECC is that there does not exist, 
at this point, a result analogous to the threshold theorem 
of fault tolerant quantum error correction. This means 
that we cannot yet guarantee full sealability of ERD as 
a stand-alone method. However, in principle it is always 
possible to concatenate ERD with QECC, as done, e.g., 
for DFS with QECC in || || II 11 > and then the 
Standard fault tolerance results apply. 



Finally, we note that ERD is a general method, that 
is not limited to trapped ions. We hope that the meth- 
ods proposed here will inspire experimentalists to im- 
plement encoded recoupling and decoupling in the lab, 
thus demonstrating the possibility of fully decoherence- 
protected quantum computation, in particular using 
trapped ions. 
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